Inherent Structures and Kauzmann Temperature of Confined Liquids 
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Calculations of the thermodynamical properties of a supercooled liquid confined in a matrix are 
performed with an inherent structure analysis. The liquid entropy is computed by means of a 
thermodynamical integration procedure. The contributions to the free energy of the liquid can 
be decoupled also in confinement in the configurational and the vibrational part. We show that 
the vibrational entropy can be calculated in the harmonic approximation as in the bulk case. The 
Kauzmann temperature of the confined system is estimated from the behavior of the configurational 
entropy. 
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I. INTRODUCTION 

It is well known that most liquids upon supercooling 
undergo a transition to an amorphous state, where me- 
chanical properties typical of a solid phase combine with 
a microscopically disordered structure 0, 0, 0, @| ■ Just 
below the melting temperature the supercooled liquids 
manifest a slowing down of dynamics. This behavior 
has been successfully interpreted in terms of the Mode 
Coupling Theory (MCT) which is able to predict the 
asymptotic properties of the density correlators upon 
decreasing temperature on approaching a temperature 
Tc [3 ■ This temperature marks a crossover from a re- 
gion where the exploration of the phase space of the sys- 
tem is determined by structural relaxations, to a region 
where it is determined by hopping processes. In the last 
few years theoretical approaches based on the analysis 
of the potential energy landscape (PEL) of the super- 
cooled liquid have driven a significant progress in the 
study of the thermodynamics of the glass transition be- 
low t c n s s hes m m m m is m. 

The phenomena related to the glass transition are yet 
not well understood in the case of liquids in confined ge- 
ometries or at contact with solid surfaces although these 
situations are very relevant for many technological and 
biological applications. It is in fact still not clear how the 
theoretical approaches developed for bulk supercooled 
liquids can be extended to describe the corresponding 
phenomenology when liquids are confined. 

While it has been shown that MCT works also for 
interpreting the dynamics of confined liquids in several 
cases ll7|, [13, [ljj, |2fJ, |2l| , no studies of dynamics in con- 
finement below Tc have so far been performed. It is 
therefore very relevant to study how the PEL and the 
thermodynamical properties below Tc are modified by 
the presence of confinement. Recently a mean field anal- 
ysis of the PEL for thin films has shown that confinc- 
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ment could affect the thermodynamical behaviour and 
the glass transition 

As proposed by different authors [2^, the behavior 
of a bulk supercooled liquid is determined by the dynam- 
ics of the system in and between the basins of the PEL. At 
low enough temperature there are two separated regimes, 
the dynamics on the short time scale can be described as 
the motion around the local minima, while the long time 
dynamics is related to the transition between different 
basins of energy. This separation of reg imes has been 
framed by Stillinger and Weber (SW) |24( in the formal- 
ism of the inherent structures (IS). According to their 
definition an IS is the configuration of local minima of 
the PEL. A basin is the set of points which maps to the 
same IS under a local energy minimization performed by 
a steepest descent procedure starting from a configura- 
tion equilibrated at a certain temperature. In the SW 
formulation under the assumption that the basins with 
the same IS energy eis have equivalent properties in the 
canonical partition function the motions between differ- 
ent basins and the vibrations inside a single basin can 
be decoupled. This formulation allows one to define and 
study a configurational entropy S con f . This quantity, 
which represents the difference between the liquid and 
the disordered solid entropies, plays a central role in un- 
derstanding the glass transition. 

In the process of cooling the configurational entropy 
decreases and eventually vanishes at a finite temperature, 
defined as the the Kauzmann temperature Tk (25|. In the 
interpretation of Adam, Gibbs and Di Marzio HH, H3, H^] 
at Tk an ideal thermodynamical transition should take 
place from the supercooled liquid to an amorphous phase 
with a single configuration 0, U, EHH HI EJ The sin- 
gular behavior of thermodynamical quantities measured 
in experiments at the conventional glass transition tem- 
perature T g would be related to the true transition oc- 
curring at Tk < T g . 

It is generally found that phase transitions in confined 
fluids are modified by confinement both from geometric 
effects and the interaction with the substrate. This is 
particularly true when large spatial correlations are ex- 
pected to take place and finite size effects could influence 



2 



the transition in a fluid confined in a restricted environ- 
ment. This would be also the case for the glass transition 
at Tk if it is interpreted in terms of a second order phase 
transition |3 ill or in the framework of the mosaic state 
scenario [32| . 

Here we consider the case of a glass forming confined 
liquid, a Lennard- Jones binary mixture (LJBM), embed- 
ded in a disordered array of soft spheres. Molecular Dy- 
namics simulations have been performed for this system 
upon cooling and a numerical test of MCT properties 
has been carried out. The mixture follows also in this 
confining environment, as in the bulk, MCT predictions 
very well 0, 0] . Nonetheless important differences due 
to confinement are found. In particular the range of va- 
lidity of the MCT predictions suffers a reduction of 60% 
with respect to the bulk. We found a crossover temper- 
ature Tc = 0.356 (in Lennard-Jones units) to be com- 
pared with the bulk value Tc = 0.435 (in Lennard-Jones 
units) .3.3] and therefore we observed a reduction of circa 
20% of Tc in going from the bulk to the confined LJBM. 

We performed in this paper an IS analysis to eval- 
uate the IS distributions, the temperature dependence 
of the configurational entropy and finally the Kauz- 
mann temperature of the confined LJBM to be compared 
with the corresponding values for the LJBM in the bulk 
phase 0,[I2. The paper develops as follows: in the next 
section we report computational details. The third sec- 
tion is devoted to the calculation of the IS for the confined 
LJBM. In the fourth section we evaluate the configura- 
tional entropy and the Kauzmann temperature. The last 
section is devoted to the conclusions. 



II. COMPUTER SIMULATION OF THE 
CONFINED LIQUID UPON SUPERCOOLING 

We studied the LJBM proposed in ref. j^S] embed- 
ded in a rigid disordered array of 16 soft spheres. The 
liquid binary mixture is composed by 800 particles of 
type A and 200 particles of type B. The parameters of 
the Lennard-Jones potential are 6aa = 1, &aa = 1, 
(bb = 0.5, 0bb = 0.88, (lab = 1.5, and uab = 0.8. 
In the following the LJ units will be used. The A and B 
particles interact with the soft spheres with a potential 
V{r) = e(a/r) 12 where e SA = 0.32, <j sa = 3, e SB = 0.22, 
a SB = 2.94. 

Molecular dynamics simulations have been performed 
in the NVT ensemble along a isochoric path at various 
temperatures upon cooling. The box length is fixed to 
L = 12.6. In previous work we already investigated the 
system in the range of temperature from T = 5.0 to T = 
0.37 and further MD simulation details are reported in 
ref. [H[l|. 

Starting from the equilibrated configura- 
tions at the following temperatures: T = 
5,2,0.8,0.6,0.55,0.5,0.475,45,0.425,0.38, we per- 
formed new simulations for each temperature in order 
to obtain a number of equally spaced configurations and 



calculate the corresponding IS for each temperature. 

The IS have been obtained by the conjugate-gradient 
minimization procedure described in the literature and 
adapted to our confined system. For each temperature 
1000 configurations have been minimized. The Hessian 
matrix has been diagonalized at each IS to calculate the 
eigenfrequencies . 



III. INHERENT STRUCTURE ANALYSIS 

In the SW formulation of IS the canonical partition 
function can be written as follows: 

Z N {T) = 

J de I sn(e I s)exp{-{e IS + f(T,e IS )]/k B T} (1) 

where f2(e/s) is the number of distinct basins with en- 
ergy eis and f(T, e/s) is the free energy of the system 
restricted to a single basin with energy e/s- The con- 
figurational entropy S con f can be defined as Scon} — 
k B ln{VL{eis))- The energies of the IS are distributed with 
a probability given by 

P(eis,T) = 

exp[~ (e/s + / (T, e IS ) - TS conf (e/s)) /k B T] 

Z N {T) K ' 

The configurational entropy can be defined also as the 
difference between the liquid entropy and the entropy of 
the disordered solid (DS) 

Sconf = Sliquid Sds (3) 

From the quenches performed at each temperature we 
calculate the distribution functions P(e/5,T) of the IS. 
These are shown in Fig. ^ together with Gaussian best 
fits. The gaussian curves appear to reproduce sufficiently 
well the distribution functions. For the lowest tempera- 
tures the curves are narrower and more peaked around 
the average value. 

We can now look at the behavior of the configurational 
entropy by considering Eq. J2J), from which we obtain 

In [P {e IS ,T)}+e IS /k B T = 

Sconf (e/ S ) /k B - / (T, e/ S ) /k B T - In [Z N {T)\ (4) 

The left hand side of this equation can be calculated from 
the distribution functions of Fig.Q]to obtain a new set of 
curves. By plotting all these curves as function of e/s we 
see that they can be superimposed by subtracting a tem- 
perature dependent term as shown in Fig.|21for T < 0.80. 
Looking at the right hand side of Eq.0|this result implies 
that the the basin free energy f(T, e/s) is almost inde- 
pendent of the IS energy. The master curve represents, 
apart from an unknown temperature dependent term, the 
configurational entropy. 
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FIG. 1: Distribution functions P (e/s,T). Lower tempera- 
tures are on the left. The functions are omitted for temper- 
atures T = 0.50 and T = 0.55 since they are almost overim- 
posed to the ones at T = 0.475 and T = 0.60 respectively. 
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FIG. 2: Entropy per particle calculated from the dis- 
tribution functions of the inherent structure energy ac- 
cording to Eq. 21 a t the different temperatures: T — 
0.38, 0.425, 0.45, 0.475, 0.50, 0.55, 0.60, 0.80, 2.0, 5.0. The un- 
known temperature dependent term is obtained by maximiz- 
ing the overlap between the curves. For T > 0.80 the curves 
deviate from the master curve. 

Since the basin free energy is approximately indepen- 
dent of eis for T < 0.8 in the confined liquid, the parti- 
tion function defined in Eq. ^ can be separated as 

Z N (T) « 

exp[-f(T, e IS )/k B T] J de Is fl (e IS ) exp (-e IS /k B T)(5) 

The confined liquid at low enough temperature can 
be assumed to be composed of an inherent structure 
subsystem in thermal equilibrium with the vibrational 
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FIG. 3: Density of state (DOS) of the confined liquid com- 
pared with the DOS of the bulk liquid. The DOS are obtained 
by quenching the systems from T = 0.38 for the confined 
mixture and from T — 0.51 for the bulk. Both systems were 
equilibrated at the same pressure before the quench, the tem- 
peratures are different due to the fact that Tc is lower for the 
confined system |34|. 



subsystem. The IS represents the long time dynam- 
ics of the system due to transitions between the differ- 
ent basins of energy, whose degeneracy is counted by 
^ (e/s) = exp(S con f /ksT). The vibrational spectrum 
is related to the oscillations close to the minimum of the 
single basin. It can be obtained by diagonalizing the 
Hessian matrix after a conjugate gradient minimization 
starting from equivalent state points [34[. 

In Fig. [3] we report the comparison of the density of 
state (DOS) of the vibrational spectrum of the confined 
and the bulk LJBM obtained with the same method. 
We observe that the confinement does not induce large 
changes in both the shape and the spectral range of the 
eigenfrequencies. For the bulk LJBM it has been also 
shown that the basin free energy can be approximated 
with the harmonic vibrational contribution We will 
come back later to this point. 



IV. CONFIGURATION AL ENTROPY AND 
KAUZMANN TEMPERATURE 

The behavior of the configurational entropy has been 
determined from Eq. 21 and shown in Fig. |2 but the cal- 
culation of the Kauzmann temperature requires the ab- 
solute value of S con f. 

A thermodynamical integration procedure allows one 
to evaluate the full entropy of the liquid including the 
temperature dependent integration constant which ap- 
pears in Eq. 0] Starting from a state reference point at 
temperature T r at the given volume V of the simulation 
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box, the entropy Stot can be computed as 

Stot (T, V) = S ref (T r: V) + £ {^P-) v dT' (6) 

where U(T) is the internal energy calculated in the sim- 
ulation along the isochoric path. As reference point we 
assume T r = 5.0. The reference entropy S re f (T r , V) is 
derived from the corresponding Sb u ik(T r , V) of the bulk 
system at the same temperature T r and volume V by 
adding the contribution of the work needed to include 
the soft spheres by keeping the volume constant. 

S r ef{T r , V) = 

S b r f f k (T r , V) + 1 [U(T r , V) - U mk {T r , V)] (7) 

Sbuik(T r ,V) is obtained as follows: 

S% l f k (T r ,V) = S%*(T r ,V) 
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Ubulk(T r , V) 
T 

T 



V pbulk 



-dV 



(8) 



where P^ c is the excess pressure of the bulk and 
Sideai(Tr,V) is the entropy of the ideal two component 
gas. 

The result is the topmost curve shown in Fig.^J Below 
the lowest investigated temperature T = 0.38 the curve 
is extrapolated by an accurate polynomial fit. 

At variance with the bulk in our case the effective den- 
sity of the confined liquid is not constant, since the free 
volume accessible to the A and B particles changes with 
the temperature due to the soft spheres interaction po- 
tential [13, HI, |34( • The calculation of S to t has been per- 
formed along an isochoric path and the internal energy 
of the confined liquid used in Eq. contains also a con- 
tribution W con fi n (T) due to the work done to change the 
effective density of the liquid inside the simulation box at 
constant volume. This contribution has to be subtracted 
to extract the entropy of the liquid from which the the 
configurational entropy can be obtained 

Su q (T,V)=S tot (T,V)-£ ^<^^-^dT' (9) 

This integral can be calculated referring to a correspond- 
ing bulk system simulated at the same pressures and tem- 
peratures of the confined mixture, by considering 
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con fin 



dW, 



con fin ^ *liq 



dT 



dV, 



liq 



dT 



(10) 



where Vn q is the effective volume of the confined liquid. 
Vii q (T) can be derived by comparison with an equivalent 
bulk at the same pressure. The thermodynamical path 
followed by the equivalent bulk and the corresponding 
densities will be published in a separate paper |34|. The 
result for Sn q is also reported in Fig. 0] We note that 
Suq of the confined system when compared to that of 
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FIG. 4: Total entropy of the confined system as obtained from 
the thermodynamical integration (full line) as function of T. 
The dashed curve is the entropy Su q obtained after subtract- 
ing the contribution due to W CO nfin, see Eq. [5] In the figure 
are also reported the vibrational entropy Shds (dotted line) 
calculated in the harmonic approximation and the quantity 
Shds ~ 3ln(T/To (broken line with filled squares) fitted with 
the polynomial 3.704 + 0.1555T - 0.0835T 2 . In the inset S Hq 
and Shds are reported on a larger scale to better show the 
crossing point corresponding to the Kauzmann temperature 
Tk- All entropies are per particle and in units of fcg. 



the bulk 0, Il2| assumes higher values for high tempera- 
tures but approaches the zero value at approximatively 
the same temperature as the bulk. Therefore the en- 
hancement of entropy due to the additional disorder in- 
duced by the presence of the soft spheres seems to become 
less marked as the temperature is decreased. 

Assuming the harmonic approximation to be valid also 
for our confined LJBM, from the eigenfrequencies ob- 
tained from the IS we can evaluate the entropy of the 
harmonic disordered solid (HDS) with the formula 



Shds = ^ t 1 _ ln 



(11) 



The result is also reported in Fig.0] In the same figure 
it is shown the quantity Shds /NUb ~ 3Zn(T/To), where 
To = 1 has been chosen as a reference temperature. It 
has a very weak T-dcpcndcncc well fitted by a quadratic 
polynomial. This shows that almost all the T dependence 
of the entropy Shds is contained in the term independent 
of the frequency distribution. 

With the assumptions done for the vibrational spec- 
trum of the basins we can identify the Shds with the 
entropy of the disordered solid Sds- In this way the con- 
figurational entropy S con f is obtained from the difference 
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conf 



Sua — Shds 



(12) 



From Fig. 01 as best evident in the blow up, we find that 
S conf = at the temperature T = 0.292 ± 0.02 which 
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can be identified as the Kauzmann temperature of the 
confined LJBM. 



V. DISCUSSION AND CONCLUSIONS 

We have shown that for a LJBM the IS analysis can be 
performed also in confinement. The absolute value of the 
entropy of the confined liquid can be obtained by ther- 
modynamical integration by means of a procedure where 
one refers to an equivalent bulk system at the same tem- 
perature and volume as the confined liquid for including 
the ideal terms. The result must be corrected for the 
work done to change the density of the confined liquid 
keeping constant the volume of the simulation cell. The 
correction is calculated by comparison with a bulk liquid 
at the same pressure as the confined one. The combi- 
nation of the IS analysis and the thermodynamical in- 
tegration technique allows to determine the Kauzmann 
temperature Tk of the system defined as the tempera- 
ture at which the configurational entropy vanishes. 

With the entropy of the disordered solid calculated in 
the harmonic approximation as in the bulk, we found 
that Tk = 0.292 for the confined system to be compared 
with Tk = 0.297 for the bulk. We observe therefore 
only a slight decrease of Tk upon confinement while a 
much more marked decrease is instead detected for the 
MCT cross-over temperature Tc. We obtained in fact 
Tc = 0.356 in confinement against Tc — 0.435 in the 

bulk dial. 

These results seem to confirm the connection between 
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